Method for Controlling Steady Flight of Unmanned Aircraft

ABSTRACT

Disclosed is a method for controlling stable flight of an unmanned aircraft, comprising the following steps: acquiring real-time flight operation data of the aircraft itself by means of an attitude sensor, a position sensor and an altitude sensor mounted to the unmanned aircraft, performing corresponding analysis on a kinematic problem of the aircraft by a processor mounted thereto, and establishing a dynamics model of the aircraft (S1); designing a controller of the unmanned aircraft according to a multi-layer zeroing neurodynamic method (S2); solving output control quantities of motors of the aircraft by the designed multi-layer zeroing neural network controller using the acquired real-time operation data of the aircraft and target attitude data (S3); and transferring solution results to a motor governor of the aircraft, and controlling powers of the motors according to a relationship between the control quantities solved by the controller and the powers of the motors of the multi-rotor unmanned aircraft, so as to control the motion of the unmanned aircraft (S4). Based on the multi-layer zeroing neurodynamic method, a correct solution to the problem can be approached rapidly, accurately and in real time, and a time-varying problem can be significantly solved.

TECHNICAL FIELD

The present invention relates to a flight control method, in particular to a method for controlling stable flight of an unmanned aircraft.

BACKGROUND ART

In recent years, the world's unmanned aircraft technology has developed rapidly. Multi-rotor aircrafts with vertical takeoff and landing, stable hovering, wireless transmission, long-range aerial photography and autonomous cruising capabilities have broad application prospects in military and civil fields. Due to the characteristics of excellent maneuverability, simple mechanical structure, easy deployment and easy maintenance, the small rotor-type aircrafts are widely used in the fields of aerial photography, power inspection, environmental monitoring, forest fire prevention, disaster inspection, anti-terrorism rescue, military reconnaissance, battlefield assessment, etc. With wide application of unmanned aircrafts, the design of stable and fast-responding controllers of the unmanned aircrafts has attracted the attention of many researchers. The conventional controllers of the unmanned aircrafts are all designed on the basis of PID closed-loop control algorithms and corresponding improved control algorithms. PID controllers and feedback closed-loop control systems, which are simple in design and have good control effects, have been widely used in the design of the controllers of the aircrafts. Although the PID controllers are easy to use, the PID controllers and the power allocation schemes obtained according to the PID controllers do not achieve the desired stability of the unmanned aircrafts.

SUMMARY OF THE INVENTION

An object of the present invention is, in order to overcome the deficiencies of the prior art, to provide a method for designing a stable flight controller and a power allocation scheme, which controls stable flight of an unmanned aircraft by acquiring real-time flight operation data of the aircraft by using a sensor, solving output control quantities of motors of the aircraft by means of a multi-layer zeroing neural network, and obtaining a corresponding power allocation scheme.

The object of the present invention is achieved by means of the following technical solution:

a method for controlling stable flight of an unmanned aircraft, comprising the steps of:

1) acquiring real-time flight operation data of the aircraft itself, analyzing a kinematic problem of the aircraft, and establishing a dynamics model of the aircraft;

2) constructing a deviation function according to the real-time flight operation data acquired in step 1) and target attitude data, and constructing neurodynamic equations based on the deviation function by using a multi-layer zeroing neurodynamic method, wherein the neurodynamic equations based on the deviation function corresponding to all parameters together constitute a controller of the unmanned aircraft, and output quantities solved from differential equations of the controller are output control quantities of motors of the aircraft; and

3) controlling powers of the motors according to a relationship between the output control quantities solved in step 2) and the powers of the motors of the multi-rotor unmanned aircraft to complete motion control over the unmanned aircraft, specifically:

according to a power allocation scheme for the unmanned aircraft, the control quantities solved by the controller have the following relationship with the powers of the motors of the multi-rotor unmanned aircraft:

U=WF

where U=[u₁ u₂ u₃ u₄]^(T) refers to the output control quantities of the unmanned aircraft, F=[F₁ . . . F_(j)]^(T) refers to the powers of the motors of the unmanned aircraft, j is the number of the motors of the multi-rotor unmanned aircraft, W is a power allocation matrix of the unmanned aircraft, and the matrix W has different forms depending on different structures and the number of rotors, and needs to be determined according to the structure thereof and the number of the rotors;

the corresponding powers F of the motors are obtained by means of matrix inversion or pseudo-inversion, i.e.:

F=W ⁻¹ U

if the matrix W is a square matrix and is reversible, W⁻¹ is obtained by means of an inverse operation, and if W is not a square matrix, W⁻¹ is solved by means of a corresponding pseudo-inverse operation; and the desired powers F of the motors are finally obtained, input voltages of the motors are controlled according to a relationship between the voltages and powers of the motors to control the rotational speeds of the motors, and the control over the powers of the motors is finally realized to complete stable flight control over the unmanned aircraft.

Further, the step of performing corresponding analysis on a kinematic problem of the aircraft by a processor mounted thereto specifically comprises:

defining a ground coordinate system E and a fuselage coordinate system B, and establishing a relationship E=RB between the ground coordinate system and the fuselage coordinate system by means of a transformation matrix R, where R may be expressed as

${{R =}\quad}{\quad\begin{bmatrix} {\cos\mspace{11mu}\theta\mspace{11mu}\cos\mspace{11mu}\psi} & {\cos\mspace{11mu}\theta\mspace{11mu}\sin\mspace{11mu}\psi} & {{- s}{in}\mspace{11mu}\theta} \\ {{\sin\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\cos\mspace{11mu}\psi} - {\cos\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\psi}} & \begin{matrix} {{\sin\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\sin\mspace{11mu}\psi} +} \\ {\cos\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\psi} \end{matrix} & {\sin\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\theta} \\ {{\cos\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\cos\mspace{11mu}\psi} + {\sin\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\psi}} & \begin{matrix} {{\cos\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\sin\mspace{11mu}\psi} -} \\ {\sin\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\psi} \end{matrix} & {\cos\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\theta} \end{bmatrix}}$

where ϕ is a roll angle, θ is a pitch angle, and ψ is a yaw angle;

ignoring the effect of an air resistance on the aircraft, the stress analysis (in the form of Newton-Euler) of the aircraft system in the fuselage coordinate system is as follows

${{\begin{bmatrix} {mI}_{3 \times 3} & 0 \\ 0 & I \end{bmatrix}\begin{bmatrix} \overset{.}{V} \\ \overset{.}{\omega} \end{bmatrix}} + \begin{bmatrix} {\omega \times m\; V} \\ {\omega \times I\;\omega} \end{bmatrix}} = \begin{bmatrix} F \\ \tau \end{bmatrix}$

where m is the total mass of the aircraft, I_(3×3) is a unit matrix, I is an inertia matrix, V is a linear velocity in the fuselage coordinate system, ω is an angular velocity in the fuselage coordinate system, F is a resultant external force, and τ is a resultant torque.

Further, the step of establishing a dynamics model of the aircraft specifically comprises:

according to the defined ground coordinate system E and fuselage coordinate system B, the relationship E=RB established between the two by means of the transformation matrix R and the stress analysis of the aircraft system in the fuselage coordinate system, obtaining dynamics equations of the multi-rotor aircraft as follows

$\left\{ \begin{matrix} {\overset{¨}{x} = \frac{u_{x}u_{1}}{m}} \\ {\overset{¨}{y} = \frac{u_{y}u_{1}}{m}} \\ {\overset{¨}{z} = {{- g} + \frac{\left( {\cos\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\theta} \right)u_{1}}{m}}} \\ {\overset{¨}{\phi} = \frac{{\left( {I_{y} - I_{z}} \right)\overset{.}{\theta}\overset{.}{\psi}} + {lu_{2}}}{I_{x}}} \\ {\overset{¨}{\theta} = \frac{{\left( {I_{z} - I_{x}} \right)\overset{¨}{\psi}\overset{¨}{\phi}} + {lu_{3}}}{I_{y}}} \\ {\overset{¨}{\psi} = \frac{{\left( {I_{x} - I_{y}} \right)\phi\overset{¨}{\theta}} + u_{4}}{I_{z}}} \end{matrix} \right.$

where l is an arm length; g is a gravitational acceleration; x, y, z are respectively position coordinates of the aircraft in the ground coordinate system; I_(x), I_(y), I_(z) are respectively rotational inertia of the aircraft in X, Y and Z axes; u_(x)=cos ϕ sin θ cos ψ+sin ϕ sin ψ; u_(y)=cos ϕ sin θ sin ψ−sin ϕ cos ψ; and u₁, u₂, u₃, u₄ are output control quantities.

Further, the step of designing a controller of the unmanned aircraft specifically comprises the steps of:

(1) designing a deviation function regarding the output control quantity u₁ from the vertical altitude z by means of the multi-layer zeroing neurodynamic method, and designing an altitude controller for the unmanned aircraft according to this deviation function;

(2) designing a deviation function regarding u_(x) and u_(y) from the horizontal positions x and y by means of the multi-layer zeroing neurodynamic method, designing a position controller for the unmanned aircraft according to this deviation function, and then inversely solving target attitude angles ϕ_(T) and θ_(T); and

(3) designing a deviation function regarding the output control quantities u₂˜u₄ from the roll angle ϕ, the pitch angle θ and the yaw angle ψ by means of the multi-layer zeroing neurodynamic method, and designing an attitude controller according to this deviation function.

Further, the step of designing a deviation function regarding the output control quantity u₁ and a corresponding altitude controller of the unmanned aircraft specifically consists in that:

for the vertical altitude z, according to the target altitude value z_(T) and the actual altitude value z in the Z axis, a deviation function may be defined as

e _(z1) =z−z _(T)  (1)

and its derivative may be obtained as follows

ė _(z1) =ż−ż _(T)  (2)

in order to converge the actual value z to the target value z_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(z1) =γe _(z1)  (3)

where γ is a constant;

equations (1) and (2) are substituted into equation (3) and collating is performed to obtain

ż−ż _(T)+γ(z−z _(T))=0  (4)

since equation (4) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(z2) =ż−ż _(T)+γ(z−z _(T))  (5)

and its derivative may be obtained as follows

ė _(z2) ={umlaut over (z)}−{umlaut over (z)} _(T)+γ(ż−ż _(T))  (6)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(z2) =−γe _(z2)  (7)

equations (5) and (6) are substituted into equation (7) and collating is performed to obtain

{umlaut over (z)}−{umlaut over (z)} _(T)+2γ(ż−ż _(T))+γ²(z−z _(T))=0  (8)

in this way, a deviation function may be defined as

E _(z) ={umlaut over (z)}−{umlaut over (z)} _(T)+2γ(ż−ż _(T))+γ²(z−z _(T))  (9)

according to the dynamics equations of the aircraft, (9) may be simplified into

E _(z) =a _(z) u ₁ +b _(z)  (10)

where

${a_{z} = \frac{\cos\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\theta}{m}},$

and b_(z)=−g−{umlaut over (z)}_(T)+2γ(ż−ż_(T))+γ² (z−z_(T)); and its derivative may be obtained as follows

Ė _(z) =a _(z) {dot over (u)} ₁ +{dot over (a)} _(z) u ₁ +{dot over (b)} _(z)  (11)

it is possible to use the multi-layer zeroing neurodynamic method to design

Ė _(z) =−γE _(z)  (12)

equations (10) and (11) are substituted into equation (12) and collating is performed to obtain

a _(z) {dot over (u)} ₁=γ(a _(z) u ₁ +b _(z))−{dot over (b)} _(z) −{dot over (a)} _(z) u ₁  (13).

Further, the step of designing a deviation function regarding u_(x) and u_(y) and a corresponding position controller for the unmanned aircraft specifically consists in that:

for the horizontal position x, according to the target value x_(T) and the actual value x in the X axis, a deviation function may be defined as

e _(x1) =x−x _(T)  (14)

and its derivative may be obtained as follows

ė _(x1) ={dot over (x)}−{dot over (x)} _(T)  (15)

in order to converge the actual value x to the target value x_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(x1) =−γe _(x1)  (16)

equations (14) and (15) are substituted into equation (16) and collating is performed to obtain

{dot over (x)}−{dot over (x)} _(T)+γ(x−x _(T))=0  (17)

since equation (17) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(x2) ={dot over (x)}−{dot over (x)} _(T)+γ(x−x _(T))  (18)

and its derivative may be obtained as follows

ė _(x2) ={umlaut over (x)}−{umlaut over (x)} _(T)+γ({dot over (x)}−{dot over (x)} _(T))  (19)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(x2) =−γe _(x2)  (20)

equations (18) and (19) are substituted into equation (20) and collating is performed to obtain

{umlaut over (x)}−{umlaut over (x)} _(T)+2γ({dot over (x)}−{dot over (x)} _(T))+γ²(x−x _(T))=0  (21)

in this way, a deviation function may be defined as

E _(x) ={umlaut over (x)}−{umlaut over (x)} _(T)+2γ({dot over (x)}−{dot over (x)} _(T))+γ²(x−x _(T))  (22)

according to the dynamics equations of the aircraft, equation (22) may be simplified into

E _(x) =a _(x) u _(x) +b _(x)  (23)

where

${a_{x} = \frac{u_{1}}{m}},$

and b_(x)=−{umlaut over (x)}_(T)+2γ({dot over (x)}−{dot over (x)}_(T))+γ²(x−x_(T)); and its derivative may be obtained as follows

Ė _(x) =a _(x) {dot over (u)} _(x) +{dot over (a)} _(x) u _(x) +b _(x)  (24)

it is possible to use the multi-layer zeroing neurodynamic method to design

Ė _(x) =−γE _(x)  (25)

equations (23) and (24) are substituted into equation (25) and collating is performed to obtain

a _(x) {dot over (u)} _(x)=−γ(a _(x) u _(x) +b _(x))−{dot over (b)} _(x) −{dot over (a)} _(x) u _(x)  (26)

for the horizontal position y, according to the target value y_(T) and the actual value y in the Y axis, a deviation function may be defined as

e _(y1) =y−y _(T)  (27)

and its derivative may be obtained as follows

ė _(y1) =y−y _(T)  (28)

in order to converge the actual value y to the target value y_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(y1) =−γe _(y1)  (29)

equations (27) and (28) are substituted into equation (29) and collating is performed to obtain

{dot over (y)}−{dot over (y)} _(T)+γ(y−y _(T))=0  (30)

since equation (30) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(y2) ={dot over (y)}−{dot over (y)} _(T)+γ(y−y _(T))  (31)

and its derivative may be obtained as follows

ė _(y2) =ÿ−ÿ _(T)+γ({dot over (y)}−{dot over (y)} _(T))  (32)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(y2) =−γe _(y2)  (33)

equations (31) and (32) are substituted into equation (33) and collating is performed to obtain

ÿ−ÿ _(T)+2γ({dot over (y)}−{dot over (y)} _(T))+γ²(y−y _(T))=0  (34)

in this way, a deviation function may be defined as

E _(y) =ÿ−ÿ _(T)+2γ({dot over (y)}−{dot over (y)} _(T))+γ²(y−y _(T))  (35)

according to the dynamics equations of the aircraft, equation (35) may be simplified into

E _(y) =a _(y) u _(y) +b _(y)  (36)

where

${a_{y} = \frac{u_{1}}{m}},$

and b_(y)=ÿ_(T)+2γ({dot over (y)}−{dot over (y)}_(T))+γ²(y−y_(T)); and its derivative may be obtained as follows

Ė _(y) =a _(y) {dot over (u)} _(y) +{dot over (a)} _(y) u _(y) +{dot over (b)} _(y)  (37)

it is possible to use the multi-layer zeroing neurodynamic method to design

Ė _(y) =−γE _(y)  (38)

equations (36) and (37) are substituted into equation (38) and collating is performed to obtain

a _(y) {dot over (u)} _(y)=γ(a _(y) u _(y) +b _(y))−{dot over (b)} _(y) −{dot over (a)} _(y) u _(y)  (39).

Further, according to the designed position controller, the calculation method of inversely solving the target attitude angles ϕ_(T) and θ_(T) is:

u_(x) and u_(y) solved from equations (26) and (39) are

$\left\{ \begin{matrix} {u_{x} = {{\cos\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\cos\mspace{11mu}\psi} + {\sin\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\psi}}} \\ {u_{y} = {{\cos\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\sin\mspace{11mu}\psi} - {\sin\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\psi}}} \end{matrix} \right.$

so that the inversely solved target angle values ϕ_(T) and θ_(T) are

$\begin{matrix} {\left\{ \begin{matrix} {\phi_{T} = {\sin^{- 1}\left( {{u_{x\mspace{11mu}}\sin\mspace{11mu}\psi} - {u_{y}\mspace{14mu}\cos\mspace{11mu}\psi}} \right)}} \\ {\theta_{T} = {\sin^{- 1}\left( \frac{u_{x} - {\sin\mspace{11mu}\phi_{T}\mspace{14mu}\sin\mspace{11mu}\psi}}{\cos\mspace{11mu}\phi_{T}\mspace{14mu}\cos\mspace{11mu}\psi} \right)}} \end{matrix} \right..} & (40) \end{matrix}$

Further, the step of designing a deviation function regarding the output control quantity u₂˜u₄ and a corresponding attitude controller of the unmanned aircraft specifically consists in that:

for the roll angle ϕ, according to the target angle ϕ_(T) solved in (40) and the actual angle ϕ, a deviation function may be defined as

e _(ϕ1)=ϕ−ϕ_(T)  (41)

and its derivative may be obtained as follows

ė _(ϕ1)={dot over (ϕ)}−{dot over (ϕ)}_(T)  (42)

in order to converge the actual value ϕ to the target value ϕ_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(ϕ1) =−γe _(ϕ1)  (43)

equations (41) and (42) are substituted into equation (43) and collating is performed to obtain

{dot over (ϕ)}−{dot over (ϕ)}_(T)+γ(ϕ−ϕ_(T))=0  (44)

since equation (44) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(ϕ2)={dot over (ϕ)}−{dot over (ϕ)}_(T)+γ(ϕ−ϕ_(T))  (45)

and its derivative may be obtained as follows

ė _(ϕ2)={umlaut over (ϕ)}−{umlaut over (ϕ)}_(T)+γ({dot over (ϕ)}−{dot over (ϕ)}_(T))  (46)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(ϕ2) =−γe _(ϕ2)  (47)

equations (45) and (46) are substituted into equation (47) and collating is performed to obtain

{umlaut over (ϕ)}−{umlaut over (ϕ)}_(T)+2γ({dot over (ϕ)}−{dot over (ϕ)}_(T))+γ²(ϕ−ϕ_(T))=0  (48)

in this way, a deviation function may be defined as

E _(ϕ)={umlaut over (ϕ)}−{umlaut over (ϕ)}_(T)+2γ({dot over (ϕ)}−{dot over (ϕ)}_(T))+γ²(ϕ−ϕ_(T))  (49)

according to the dynamics equations of the aircraft, equation (49) may be simplified into

E _(ϕ) =a _(ϕ) u ₂ +b _(ϕ)  (50)

where

${a_{\phi} = \frac{l}{I_{x}}},{{{{and}\mspace{14mu} b_{\phi}} = {\frac{\left( {I_{y} - I_{z}} \right)\overset{.}{\theta}\overset{.}{\psi}}{I_{x}} - {\overset{¨}{\phi}}_{T} + {2{\gamma\left( {\overset{.}{\phi} - {\overset{.}{\phi}}_{T}} \right)}} + {\gamma^{2}\left( {\phi - \phi_{T}} \right)}}};}$

and its derivative may be obtained as follows

Ė _(ϕ) =a _(ϕ) {dot over (u)} ₂ +{dot over (a)} _(ϕ) u ₂ +{dot over (b)} _(ϕ)  (51)

according to the multi-layer zeroing neurodynamic method, it is possible to design

Ė _(ϕ) =γE _(ϕ)  (52)

equations (50) and (51) are substituted into equation (52) and collating is performed to obtain

a _(ϕ) {dot over (u)} ₂=−γ(a _(ϕ) u ₂ +b _(ϕ))−{dot over (b)} _(ϕ) −{dot over (a)} _(ϕ) u ₂  (53)

for the pitch angle θ, according to the target angle θ_(T) solved in (40) and the actual angle θ, a deviation function may be defined as

e _(θ1)=θ−θ_(T)  (54)

and its derivative may be obtained as follows

ė _(θ1)={dot over (θ)}−{dot over (θ)}_(r)  (55)

in order to converge the actual value θ to the target value θ_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(θ1) =−γe _(θ1)  (56)

equations (54) and (55) are substituted into equation (56) and collating is performed to obtain

{dot over (θ)}−{dot over (θ)}_(T)+γ(θ−θ_(T))=0  (57)

since equation (57) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(θ2)={dot over (θ)}−{dot over (θ)}_(T)+γ(θ−θ_(T))  (58)

and its derivative may be obtained as follows

ė _(θ2)={umlaut over (θ)}−{umlaut over (θ)}_(T)+γ({dot over (θ)}−{dot over (θ)}_(T))  (59)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(θ2) =−γe _(θ2)  (60)

equations (58) and (59) are substituted into equation (60) and collating is performed to obtain

{umlaut over (θ)}−{umlaut over (θ)}_(T)+2γ({dot over (θ)}−{dot over (θ)}_(T))+γ²(θ−θ_(T))=0  (61)

in this way, a deviation function may be defined as

E _(θ)={umlaut over (θ)}−{umlaut over (θ)}_(T)+2γ({dot over (θ)}−{dot over (θ)}_(T))+γ²(θ−θ_(T))  (62)

according to the dynamics equations of the aircraft, equation (62) may be simplified into

E _(θ) =a _(θ) u ₃ +b _(θ)  (63)

where

${a_{\theta} = \frac{l}{I_{y}}},{{{{and}\mspace{14mu} b_{\theta}} = {\frac{\left( {I_{z} - I_{x}} \right)\overset{.}{\phi}\overset{.}{\psi}}{I_{y}} - {\overset{¨}{\theta}}_{T} + {2{\gamma\left( {\overset{.}{\theta} - {\overset{.}{\theta}}_{T}} \right)}} + {\gamma^{2}\left( {\theta - \theta_{T}} \right)}}};}$

and its derivative may be obtained as follows

Ė _(θ) =a _(θ) {dot over (u)} ₃ +{dot over (b)} _(θ)  (64)

according to the multi-layer zeroing neurodynamic method, it is possible to design

Ė _(θ) =−γE _(θ)  (65)

equations (63) and (64) are substituted into equation (65) and collating is performed to obtain

a _(θ) {dot over (u)} ₃=γ(a _(θ) u ₃ +b _(θ))−{dot over (b)} _(θ) −{dot over (a)} _(θ) u ₃  (66)

for the yaw angle ψ, according to an artificially set angle ψ_(T) and the actual angle ψ, a deviation function may be defined as

e _(ψ1)=ψ−ψ_(T)  (67)

and its derivative may be obtained as follows

ė _(ψ1)={dot over (ψ)}−{dot over (ψ)}_(T)  (68)

in order to converge the actual value ψ to the target value ψ_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(ψ1) =−γe _(ψ1)  (69)

equations (67) and (68) are substituted into equation (69) and collating is performed to obtain

{dot over (ψ)}−{dot over (ψ)}_(T)+γ(ψ−ψ_(T))=0  (70)

since equation (70) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(ψ2)={dot over (ψ)}−{dot over (ψ)}_(T)+γ(ψ−ψ_(T))  (71)

and its derivative may be obtained as follows

ė _(ψ2)={umlaut over (ψ)}−{umlaut over (ψ)}_(T)+γ({dot over (ψ)}−{dot over (ψ)}_(T))  (72)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(ψ2) =−γe _(ψ2)  (73)

equations (71) and (72) are substituted into equation (73) and collating is performed to obtain

{umlaut over (ψ)}−{umlaut over (ψ)}_(T)+2γ({dot over (ψ)}−{dot over (ψ)}_(T))+γ²(ψ−ψ_(T))=0  (74)

in this way, a deviation function may be defined as

E _(ψ)={umlaut over (ψ)}−{umlaut over (ψ)}_(T)+2γ({dot over (ψ)}−{dot over (ψ)}_(T))+γ²(ψ−ψ_(T))  (75)

according to the dynamics equations of the aircraft, equation (75) may be simplified into

E _(ψ) =a _(ψ) u ₄ +b _(ψ)  (76)

where

${a_{\psi} = \frac{1}{I_{z}}},{{{{and}\mspace{14mu} b_{\psi}} = {\frac{\left( {I_{x} - I_{y}} \right)\overset{.}{\phi}\overset{.}{\theta}}{I_{z}} - {\overset{¨}{\psi}}_{T} + {2{\gamma\left( {\overset{.}{\psi} - {\overset{.}{\psi}}_{T}} \right)}} + {\gamma^{2}\left( {\psi - \psi_{T}} \right)}}};}$

and its derivative may be obtained as follows

Ė _(ψ) =a _(ψ) {dot over (u)} ₄ +{dot over (a)} _(ψ) u ₄ +{dot over (b)} _(ψ)  (77)

according to the multi-layer zeroing neurodynamic method, it is possible to design

Ė _(ψ) =−γE _(ψ)  (78)

equations (76) and (77) are substituted into equation (78) and collating is performed to obtain

a _(ψ) {dot over (u)} ₄=−γ(a _(ψ) u ₄ +b _(ψ))−{dot over (b)} _(ψ) −{dot over (a)} _(ψ) u ₄  (79).

Further, the step in which the altitude controller designed according to the altitude variable z, the position controller designed according to the position variables x and y, and the attitude controller designed according to the attitude control quantities ϕ, θ and ψ together constitute a stable flight controller of the multi-rotor unmanned aircraft specifically consists in that:

a controller of the unmanned aircraft may be obtained according to equations (13), (53), (66) and (79), wherein the controller can be implemented by a network structure; the controller of the unmanned aircraft is capable of controlling the stable flight of the unmanned aircraft; and the controller may be written in the following form:

$\quad\left\{ \begin{matrix} {{\overset{.}{u}}_{1} = \frac{{- {\gamma\left( {{a_{z}u_{1}} + b_{z}} \right)}} - {\overset{.}{b}}_{z} - {{\overset{.}{a}}_{z}u_{1}}}{a_{z}}} \\ {{\overset{.}{u}}_{2} = \frac{{- {\gamma\left( {{a_{\phi}u_{2}} + b_{\phi}} \right)}} - {\overset{.}{b}}_{\phi} - {{\overset{.}{a}}_{\phi}u_{2}}}{a_{\phi}}} \\ {{\overset{.}{u}}_{3} = \frac{{- {\gamma\left( {{a_{\theta}u_{3}} + b_{\theta}} \right)}} - {\overset{.}{b}}_{\theta} - {{\overset{.}{a}}_{\theta}u_{3}}}{a_{\theta}}} \\ {{\overset{.}{u}}_{4} = \frac{{- {\gamma\left( {{a_{\psi}u_{4}} + b_{\psi}} \right)}} - {\overset{.}{b}}_{\psi} - {{\overset{.}{a}}_{\psi}u_{4}}}{a_{\psi}}} \end{matrix} \right.$

a zeroing neural network is constructed from the differential equations of the controller, and the control quantities of the unmanned aircraft are solved by means of the zeroing neural network.

Compared with the prior art, the present invention has the following beneficial effects:

1. The multi-layer zeroing neural network has better convergence characteristics, can realize real-time response of the aircraft and has a strong robustness, and the controller system designed according to the neural network is stable and has a good control effect.

2. The present invention is based on the multi-layer zeroing neurodynamic method, the method is described by using a ubiquitous implicit dynamics model, derivative information of various time-varying parameters can be fully utilized from the method and system level, and the method has a certain predictive ability for solving problems, can quickly, accurately and timely approach correct solutions of the problems, and can solve a variety of time-varying problems such as matrices, vectors, algebras and optimization.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a flow chart of a method for controlling stable flight of a multi-rotor aircraft according to an embodiment of the present invention;

FIG. 2 is a side view showing the structure of the multi-rotor aircraft according to the present invention;

FIG. 3 is a top view showing the structure of the multi-rotor aircraft according to the present invention;

FIG. 4 is a three-dimensional view showing the structure of the multi-rotor aircraft according to the present invention; and

FIG. 5 is a diagram showing a fuselage coordinate system of the multi-rotor aircraft.

DETAILED DESCRIPTION OF EMBODIMENTS

Hereafter the present invention will be further described in detail in conjunction with embodiments and appended drawings, but the embodiments of the present invention are not limited thereto.

Embodiment

As shown in FIG. 1, the present embodiment provides a method for controlling stable flight of an unmanned aircraft. The method comprises the following steps:

S1: acquiring real-time flight operation data of the aircraft itself by means of an attitude sensor, a position sensor and an altitude sensor mounted to the unmanned aircraft, performing corresponding analysis on a kinematic problem of the aircraft by a processor mounted thereto, and establishing a dynamics model of the aircraft;

One type of rotor flight structure in the multi-rotor aircraft is shown in FIGS. 2, 3 and 4. The structure is a six-rotor aircraft mechanism model consisting of multi-rotor aircraft propellers, brushless motors, rotor arms and a fuselage. Arrows in FIGS. 3 and 4 indicate the directions of rotation of the motors, and the combination of the illustrated clockwise and counterclockwise directions of rotation is to achieve mutual offsetting of torques the motors so as to achieve stable steering control.

Real-time attitude data θ(t), ϕ(t) and ψ(t) of the aircraft may be acquired by sensors such as gyros and accelerometers mounted to the multi-rotor aircraft by means of quaternion algebra, Kalman filtering and other algorithms, and position data x(t), y(t) and z(t) of the aircraft in the three-dimensional space is acquired by using altitude sensors and position sensors.

The definition of aircraft attitude variables is shown in FIG. 5.

The multi-rotor aircraft in FIG. 5 is defined as follows based on the fuselage coordinate system:

(1) six motors of the six-rotor aircraft are defined No. 1 to No. 6 in the clockwise direction;

(2) X axis extends in the direction of No. 1 rotor arm and points to the forward direction of the aircraft through the center of gravity of the fuselage;

(3) Y axis extends in the direction of the axis of symmetry of No. 2 and No. 3 rotor arms and points to the right motion direction of the aircraft through the center of gravity of the fuselage;

(4) Z axis extends upwardly perpendicular to the plane of six rotors and points to the climbing direction of the aircraft through the center of gravity of the fuselage;

(5) the pitch angle θ is an angle between the X axis of the fuselage and the horizontal plane, and is set to be positive when the fuselage is downward;

(6) the roll angle φ is an angle between the Z axis of the fuselage and the vertical plane passing through the X axis of the fuselage, and is set to be positive when the fuselage is rightward; and

(7) the yaw angle ψ is an angle between the projection of the X axis of the fuselage on the horizontal plane and the X axis of a geodetic coordinate system, and is set to be positive when the head of the aircraft is leftward.

According to different rotor-type aircraft models, physical model equations and dynamics equations for the aircraft are established, and dynamics analysis may be completed by means of the following aircraft dynamics modeling steps:

defining a ground coordinate system E and a fuselage coordinate system B, and establishing a relationship E=RB between the ground coordinate system and the fuselage coordinate system by means of a transformation matrix R, where R may be expressed as

$R = \begin{bmatrix} {\cos\;{\theta cos}\psi} & {\cos\;{\theta sin}\psi} & {{- \sin}\;\theta} \\ {{\sin\;{\phi sin}\theta\cos\;\psi} - {\cos\;{\phi sin}\psi}} & {{\sin\;{\phi sin}\theta\sin\;\psi} + {\cos\;{\phi cos}\psi}} & {\sin\;{\phi cos}\theta} \\ {{\cos\;{\phi sin}\theta\cos\;\psi} + {\sin\;{\phi sin}\psi}} & {{\cos\;{\phi sin}\theta\sin\;\psi} - {\sin\;{\phi cos}\psi}} & {\cos\;{\phi cos}\theta} \end{bmatrix}$

where ϕ is a roll angle, θ is a pitch angle, and ψ is a yaw angle;

ignoring the effect of an air resistance on the aircraft, the stress analysis (in the form of Newton-Euler) of the aircraft system in the fuselage coordinate system is as follows

${{\begin{bmatrix} {mI}_{3 \times 3} & 0 \\ 0 & I \end{bmatrix}\begin{bmatrix} \overset{.}{V} \\ \overset{.}{\omega} \end{bmatrix}} + \begin{bmatrix} {\omega \times mV} \\ {\omega \times I\;\omega} \end{bmatrix}} = \begin{bmatrix} F \\ \tau \end{bmatrix}$

where m is the total mass of the aircraft, I_(3×3) is a unit matrix, I is an inertia matrix, V is a linear velocity in the fuselage coordinate system, ω is an angular velocity in the fuselage coordinate system, F is a resultant external force, and τ is a resultant torque.

According to the above equation, the dynamics equations of the aircraft can be obtained as follows

$\quad\left\{ \begin{matrix} {\overset{¨}{x} = \frac{u_{x}u_{1}}{m}} \\ {\overset{¨}{y} = \frac{u_{y}u_{1}}{m}} \\ {\overset{¨}{z} = {{- g} + \frac{\left( {\cos\phi\cos\theta} \right)u_{1}}{m}}} \\ {\overset{¨}{\phi} = \frac{{\left( {I_{y} - I_{z}} \right)\overset{.}{\theta}\overset{.}{\psi}} + {lu_{2}}}{I_{x}}} \\ {\overset{¨}{\theta} = \frac{{\left( {I_{z} - I_{x}} \right)\overset{¨}{\psi}\overset{¨}{\phi}} + {lu_{3}}}{I_{y}}} \\ {\overset{¨}{\psi} = \frac{{\left( {I_{x} - I_{y}} \right)\overset{¨}{\phi}\overset{¨}{\theta}} + u_{4}}{I_{z}}} \end{matrix} \right.$

where l is an arm length; g is a gravitational acceleration; x, y, z are respectively position coordinates of the aircraft in the ground coordinate system; I_(x), I_(y), I_(z) are respectively rotational inertia of the aircraft in X, Y and Z axes; u_(x)=cos ϕ sin θ cos ψ+sin ϕ sin ψ; u_(y)=cos ϕ sin θ sin ψ−sin ϕ cos ψ; and u₁, u₂, u₃, u₄ are output control quantities.

S2: designing a controller of the unmanned aircraft according to a multi-layer zeroing neurodynamic method;

A deviation function regarding the output control quantity u₁ is designed from the vertical altitude z, an altitude controller of the multi-rotor unmanned aircraft is designed according to this deviation function, u₁ is solved; a deviation function regarding u_(x) and u_(y) and a corresponding position controller for the multi-rotor unmanned aircraft are designed from the horizontal positions x and y, and target attitude angles ϕ_(T) and θ_(T) are inversely solved; and a deviation function regarding the output control quantities u₂˜u₄ is designed from the roll angle ϕ, the pitch angle θ and the yaw angle ψ according to the target attitude angles, and a corresponding multi-layer zeroing neural network controller is designed. The specific steps are as follows:

for the vertical altitude z, according to the target altitude value z_(T) and the actual altitude value z in the Z axis, a deviation function may be defined as

e _(z1) =z−z _(T)  (1)

and its derivative may be obtained as follows

ė _(z1) =ż−ż _(T)  (2)

in order to converge the actual value z to the target value z_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(z1) =−γe _(z1)  (3)

where γ is a constant;

equations (1) and (2) are substituted into equation (3) and collating is performed to obtain

ż−ż _(T)+γ(z−z _(T))=0  (4)

since equation (4) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(z2) =ż−ż _(T)+γ(z−z _(T))  (5)

and its derivative may be obtained as follows

ė _(z2) ={umlaut over (z)}−{umlaut over (z)} _(T)+γ(ż−ż _(T))  (6)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(z2) =−γe _(z2)  (7)

equations (5) and (6) are substituted into equation (7) and collating is performed to obtain

{umlaut over (z)}−{umlaut over (z)} _(T)+2γ(ż−ż _(T))+γ²(z−z _(T))=0  (8)

in this way, a deviation function may be defined as

E _(z) ={umlaut over (z)}−{umlaut over (z)} _(T)+2γ(ż−ż _(T))+γ²(z−z _(T))  (9)

according to the dynamics equations of the aircraft, (9) may be simplified into

E _(z) =a _(z) u ₁ +b _(z)  (10)

where

${a_{z} = \frac{\cos\phi\cos\theta}{m}},$

and b_(z)=−g−{umlaut over (z)}_(T)+2γ(ż−ż_(T))+γ² (z−z_(T)); and its derivative may be obtained as follows

Ė _(z) =a _(z) {dot over (u)} ₁ +{dot over (a)} _(z) u ₁ +{dot over (b)} _(z)  (11)

it is possible to use the multi-layer zeroing neurodynamic method to design

Ė _(z) =−γE _(z)  (12)

equations (10) and (11) are substituted into equation (12) and collating is performed to obtain

a _(z) {dot over (u)} ₁=−γ(a _(z) u ₁ +b _(z))−{dot over (b)} _(z) −{dot over (a)} _(z) u ₁  (13).

for the horizontal position x, according to the target value x_(T) and the actual value x in the X axis, a deviation function may be defined as

e _(x1) =x−x _(T)  (14)

and its derivative may be obtained as follows

ė _(x1) ={dot over (x)}−{dot over (x)} _(T)  (15)

in order to converge the actual value x to the target value x_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(x1) =−γe _(x1)  (16)

equations (14) and (15) are substituted into equation (16) and collating is performed to obtain

{dot over (x)}−{dot over (x)} _(T)+γ(x−x _(T))=0  (17)

since equation (17) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(x2) ={dot over (x)}−{dot over (x)} _(T)+γ(x−x _(T))  (18)

and its derivative may be obtained as follows

ė _(x2) ={umlaut over (x)}−{umlaut over (x)} _(T)+γ({dot over (x)}−{dot over (x)} _(T))  (19)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(x2) =−γe _(x2)  (20)

equations (18) and (19) are substituted into equation (20) and collating is performed to obtain

{umlaut over (x)}−{umlaut over (x)} _(T)+2γ({dot over (x)}−{dot over (x)} _(T))+γ²(x−x _(T))=0  (21)

in this way, a deviation function may be defined as

E _(x) ={umlaut over (x)}−{umlaut over (x)} _(T)+2γ({dot over (x)}−{dot over (x)} _(T))+γ²(x−x _(T))  (22)

according to the dynamics equations of the aircraft, equation (22) may be simplified into

E _(x) =a _(x) u _(x) +b _(x)  (23)

where

${a_{x} = \frac{u_{1}}{m}},$

and b_(x)=−{umlaut over (x)}_(T)+2γ({dot over (x)}−{dot over (x)}_(T))+γ² (x−x_(T)); and its derivative may be obtained as follows

Ė _(x) =a _(x) {dot over (u)} _(x) +{dot over (a)} _(x) u _(x) +b _(x)  (24)

it is possible to use the multi-layer zeroing neurodynamic method to design

Ė _(x) =−γE _(x)  (25)

equations (23) and (24) are substituted into equation (25) and collating is performed to obtain

a _(x) {dot over (u)} _(x)=−γ(a _(x) u _(x) +b _(x))−{dot over (b)} _(x) −{dot over (a)} _(x) u _(x)  (26)

for the horizontal position y, according to the target value y_(T) and the actual value y in the Y axis, a deviation function may be defined as

e _(y1) =y−y _(T)  (27)

and its derivative may be obtained as follows

ė _(y1) =y−y _(T)  (28)

in order to converge the actual value y to the target value y_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(y1) =−γe _(y1)  (29)

equations (27) and (28) are substituted into equation (29) and collating is performed to obtain

{dot over (y)}−{dot over (y)} _(T)+γ(y−y _(T))=0  (30)

since equation (30) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(y2) ={dot over (y)}−{dot over (y)} _(T)+γ(y−y _(T))  (31)

and its derivative may be obtained as follows

ė _(y2) =ÿ−ÿ _(T)+γ({dot over (y)}−{dot over (y)} _(T))  (32)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(y2) =−γe _(y2)  (33)

equations (31) and (32) are substituted into equation (33) and collating is performed to obtain

ÿ−ÿ _(T)+2γ({dot over (y)}−{dot over (y)} _(T))+γ²(y−y _(T))=0  (34)

in this way, a deviation function may be defined as

E _(y) =ÿ−ÿ _(T)+2γ({dot over (y)}−{dot over (y)} _(T))+γ²(y−y _(T))  (35)

according to the dynamics equations of the aircraft, equation (35) may be simplified into

E _(y) =a _(y) u _(y) +b _(y)  (36)

where

${a_{y} = \frac{u_{1}}{m}},$

b_(y)=−ÿ_(T)+2γ({dot over (y)}−{dot over (y)}_(T))+γ² (y−y_(T)); and its derivative may be obtained as follows

Ė _(y) =a _(y) {dot over (u)} _(y) +{dot over (a)} _(y) u _(y) +{dot over (b)} _(y)  (37)

it is possible to use the multi-layer zeroing neurodynamic method to design

Ė _(y) =−γE _(y)  (38)

equations (36) and (37) are substituted into equation (38) and collating is performed to obtain

a _(y) {dot over (u)} _(y)=γ(a _(y) u _(y) +b _(y))−{dot over (b)} _(y) −{dot over (a)} _(y) u _(y)  (39)

u_(x) and u_(y) may be solved from equations (26) and (39),

$\quad\left\{ \begin{matrix} {u_{x} = {{\cos\phi\sin\theta\cos\psi} + {\sin\phi\sin\psi}}} \\ {u_{y} = {{\cos\phi\sin\theta\sin\psi} - {\sin\phi\cos\psi}}} \end{matrix} \right.$

so that the inversely solved target angle values ϕ_(T) and θ_(T) may be

$\begin{matrix} \left\{ \begin{matrix} {\phi_{T} = {\sin^{- 1}\left( {{u_{x}\sin\psi} - {u_{y}\cos\psi}} \right)}} \\ {\theta_{T} = {\sin^{- 1}\left( \frac{u_{x} - {\sin\phi_{T}\sin\psi}}{\cos\phi_{T}\cos\psi} \right)}} \end{matrix} \right. & (40) \end{matrix}$

for the roll angle ϕ, according to the target angle ϕ_(T) solved in (40) and the actual angle ϕ, a deviation function may be defined as

e _(ϕ1)=ϕ−ϕ_(T)  (41)

and its derivative may be obtained as follows

ė _(ϕ1)={dot over (ϕ)}−{dot over (ϕ)}_(T)  (42)

in order to converge the actual value ϕ to the target value ϕ_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(ϕ1) =−γe _(ϕ1)  (43)

equations (41) and (42) are substituted into equation (43) and collating is performed to obtain

{dot over (ϕ)}−{dot over (ϕ)}_(T)+γ(ϕ−ϕ_(T))=0  (44)

since equation (44) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(ϕ2)={dot over (ϕ)}−{dot over (ϕ)}_(T)+γ(ϕ−ϕ_(T))  (45)

and its derivative may be obtained as follows

ė _(ϕ2)={umlaut over (ϕ)}−{umlaut over (ϕ)}_(T)+γ({dot over (ϕ)}−{dot over (ϕ)}_(T))  (46)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(ϕ2) =−γe _(ϕ2)  (47)

equations (45) and (46) are substituted into equation (47) and collating is performed to obtain

{umlaut over (ϕ)}−{umlaut over (ϕ)}_(T)+2γ({dot over (ϕ)}−{dot over (ϕ)}_(T))+γ²(ϕ−ϕ_(T))=0  (48)

in this way, a deviation function may be defined as

E _(ϕ)={umlaut over (ϕ)}−{umlaut over (ϕ)}_(T)+2γ({dot over (ϕ)}−{dot over (ϕ)}_(T))+γ²(ϕ−ϕ_(T))  (49)

according to the dynamics equations of the aircraft, equation (49) may be simplified into

E _(ϕ) =a _(ϕ) u ₂ +b _(ϕ)  (50)

where

${a_{\phi} = \frac{l}{I_{x}}},{{{{and}\mspace{14mu} b_{\phi}} = {\frac{\left( {I_{y} - I_{z}} \right)\overset{.}{\theta}\overset{.}{\psi}}{I_{x}} - {\overset{¨}{\phi}}_{T} + {2{\gamma\left( {\overset{.}{\phi} - {\overset{.}{\phi}}_{T}} \right)}} + {\gamma^{2}\left( {\phi - \phi_{T}} \right)}}};}$

and its derivative may be obtained as follows

Ė _(ϕ) =a _(ϕ) {dot over (u)} ₂ +{dot over (a)} _(ϕ) u ₂ +{dot over (b)} _(ϕ)  (51)

according to the multi-layer zeroing neurodynamic method, it is possible to design

Ė _(ϕ) =γE _(ϕ)  (52)

equations (50) and (51) are substituted into equation (52) and collating is performed to obtain

a _(ϕ) {dot over (u)} ₂=γ(a _(ϕ) u ₂ +b _(ϕ))−{dot over (b)} _(ϕ) −{dot over (a)} _(ϕ) u ₂  (53)

for the pitch angle θ, according to the target angle θ_(T) solved in (40) and the actual angle θ, a deviation function may be defined as

e _(θ1)=θ−θ_(T)  (54)

and its derivative may be obtained as follows

ė _(θ1)={dot over (θ)}−{dot over (θ)}_(r)  (55)

in order to converge the actual value θ to the target value θ_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(θ1) =−γe _(θ1)  (56)

equations (54) and (55) are substituted into equation (56) and collating is performed to obtain

{dot over (θ)}−{dot over (θ)}_(T)+γ(θ−θ_(T))=0  (57)

since equation (57) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(θ2)={dot over (θ)}−{dot over (θ)}_(T)+γ(θ−θ_(T))  (58)

and its derivative may be obtained as follows

ė _(θ2)={umlaut over (θ)}−{umlaut over (θ)}_(T)+γ({dot over (θ)}−{dot over (θ)}_(T))  (59)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(θ2) =−γe _(θ2)  (60)

equations (58) and (59) are substituted into equation (60) and collating is performed to obtain

{umlaut over (θ)}−{umlaut over (θ)}_(T)+2γ({dot over (θ)}−{dot over (θ)}_(T))+γ²(θ−θ_(T))=0  (61)

in this way, a deviation function may be defined as

E _(θ)={umlaut over (θ)}−{umlaut over (θ)}_(T)+2γ({dot over (θ)}−{dot over (θ)}_(T))+γ²(θ−θ_(T))  (62)

according to the dynamics equations of the aircraft, equation (62) may be simplified into

E _(θ) =a _(θ) u ₃ +b _(θ)  (63)

where

${a_{\theta} = \frac{l}{I_{y}}},{{{{and}\mspace{14mu} b_{\theta}} = {\frac{\left( {I_{z} - I_{x}} \right)\overset{.}{\phi}\overset{.}{\psi}}{I_{y}} - {\overset{¨}{\theta}}_{T} + {2{\gamma\left( {\overset{.}{\theta} - {\overset{.}{\theta}}_{T}} \right)}} + {\gamma^{2}\left( {\theta - \theta_{T}} \right)}}};}$

and its derivative may be obtained as follows

Ė _(θ) =a _(θ) {dot over (u)} ₃ +{dot over (b)} _(θ)  (64)

according to the multi-layer zeroing neurodynamic method, it is possible to design

Ė _(θ) =−γE _(θ)  (65)

equations (63) and (64) are substituted into equation (65) and collating is performed to obtain

a _(θ) {dot over (u)} ₃=−γ(a _(θ) u ₃ +b _(θ))−{dot over (b)} _(θ) −{dot over (a)} _(θ) u ₃  (66)

for the yaw angle ψ, according to the target angle ψ_(T) solved in (40) and the actual angle ψ, a deviation function may be defined as

e _(ψ1)=ψ−ψ_(T)  (67)

and its derivative may be obtained as follows

ė _(ψ1)={dot over (ψ)}−{dot over (ψ)}_(T)  (68)

in order to converge the actual value ψ to the target value ψ_(T), according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(ψ1) =−γe _(ψ1)  (69)

equations (67) and (68) are substituted into equation (69) and collating is performed to obtain

{dot over (ψ)}−{dot over (ψ)}_(T)+γ(ψ−ψ_(T))=0  (70)

since equation (70) is generally not established in the initial situation and does not contain relevant information of the output control quantities, and the control quantities cannot be solved, a further design is needed, and a definition is then made

e _(ψ2)={dot over (ψ)}−{dot over (ψ)}_(T)+γ(ψ−ψ_(T))  (71)

and its derivative may be obtained as follows

ė _(ψ2)={umlaut over (ψ)}−{umlaut over (ψ)}_(T)+γ({dot over (ψ)}−{dot over (ψ)}_(T))  (72)

according to the multi-layer zeroing neurodynamic method, a neurodynamic equation based on the deviation function may be designed as

ė _(ψ2) =−γe _(ψ2)  (73)

equations (71) and (72) are substituted into equation (73) and collating is performed to obtain

{umlaut over (ψ)}−{umlaut over (ψ)}_(T)+2γ({dot over (ψ)}−{dot over (ψ)}_(T))+γ²(ψ−ψ_(T))=0  (74)

in this way, a deviation function may be defined as

E _(ψ)={umlaut over (ψ)}−{umlaut over (ψ)}_(T)+2γ({dot over (ψ)}−{dot over (ψ)}_(T))+γ²(ψ−ψ_(T))  (75)

according to the dynamics equations of the aircraft, equation (75) may be simplified into

E _(ψ) =a _(ψ) u ₄ +b _(ψ)  (76)

where

${a_{\psi} = \frac{1}{I_{z}}},{{{{and}\mspace{14mu} b_{\psi}} = {\frac{\left( {I_{x} - I_{y}} \right)\overset{.}{\phi}\overset{.}{\theta}}{I_{z}} - {\overset{¨}{\psi}}_{T} + {2{\gamma\left( {\overset{.}{\psi} - {\overset{.}{\psi}}_{T}} \right)}} + {\gamma^{2}\left( {\psi - \psi_{T}} \right)}}};}$

and its derivative may be obtained as follows

Ė _(ψ) =a _(ψ) {dot over (u)} ₄ +{dot over (a)} _(ψ) u ₄ +{dot over (b)} _(ψ)  (77)

according to the multi-layer zeroing neurodynamic method, it is possible to design

Ė _(ψ) =−γE _(ψ)  (78)

equations (76) and (77) are substituted into equation (78) and collating is performed to obtain

a _(ψ) {dot over (u)} ₄=−γ(a _(ψ) u ₄ +b _(ψ))−{dot over (b)} _(ψ) −{dot over (a)} _(ψ) u ₄  (79).

S3: solving output control quantities of motors of the aircraft by the designed multi-layer zeroing neural network controller using the acquired real-time operation data of the aircraft and target attitude data; and

A controller of the unmanned aircraft may be obtained according to multi-layer zeroing neural network equations (13), (53), (66) and (79), wherein the controller can be implemented by a network structure; the controller of the unmanned aircraft is capable of controlling the stable flight of the unmanned aircraft; and the controller may be written in the following form:

$\quad\left\{ \begin{matrix} {{\overset{.}{u}}_{1} = \frac{{- {\gamma\left( {{a_{z}u_{1}} + b_{z}} \right)}} - {\overset{.}{b}}_{z} - {{\overset{.}{a}}_{z}u_{1}}}{a_{z}}} \\ {{\overset{.}{u}}_{2} = \frac{{- {\gamma\left( {{a_{\phi}u_{2}} + b_{\phi}} \right)}} - {\overset{.}{b}}_{\phi} - {{\overset{.}{a}}_{\phi}u_{2}}}{a_{\phi}}} \\ {{\overset{.}{u}}_{3} = \frac{{- {\gamma\left( {{a_{\theta}u_{3}} + b_{\theta}} \right)}} - {\overset{.}{b}}_{\theta} - {{\overset{.}{a}}_{\theta}u_{3}}}{a_{\theta}}} \\ {{\overset{.}{u}}_{4} = \frac{{- {\gamma\left( {{a_{\psi}u_{4}} + b_{\psi}} \right)}} - {\overset{.}{b}}_{\psi} - {{\overset{.}{a}}_{\psi}u_{4}}}{a_{\psi}}} \end{matrix} \right.$

a zeroing neural network is constructed from the differential equations of the controller, and the control quantities of the unmanned aircraft are solved by means of the zeroing neural network.

S4: transferring solution results of step S3 to a motor governor of the aircraft, and controlling powers of the motors according to a relationship between the control quantities solved by the controller and the powers of the motors of the multi-rotor unmanned aircraft, so as to control the motion of the unmanned aircraft;

according to a power allocation scheme for the unmanned aircraft, the control quantities solved by the controller have the following relationship with the powers of the motors of the multi-rotor unmanned aircraft:

U=WF

where U=[u₁ u₂ u₃ u₄]T refers to the control quantities of the unmanned aircraft, F=[F₁ . . . F_(j)]^(T) refers to the powers of the motors of the unmanned aircraft, j is the number of the motors of the multi-rotor unmanned aircraft, and W is a power allocation matrix of the unmanned aircraft.

In order to obtain the power required by the corresponding motor, the corresponding powers of the motors F may be obtained by means of matrix inversion or pseudo-inversion, i.e.

F=W ⁻¹ U

if the matrix W is a square matrix and is reversible, W⁻¹ is obtained by means of an inverse operation, and if W is not a square matrix, W⁻¹ is solved by means of a corresponding pseudo-inverse operation; and the desired powers F of the motors are finally obtained, input voltages of the motors are controlled according to a relationship between the voltages and powers of the motors to control the rotational speeds of the motors, and the control over the powers of the motors is finally realized to complete stable flight control over the unmanned aircraft. Since different numbers and structures of the rotors affect the control mode of the multi-rotor unmanned aircraft, the matrix W has different forms depending on the structure and the number of the rotors.

Taking the six-rotor unmanned aircraft as an example, the power allocation thereof has the following relationship:

$\left\{ \begin{matrix} {u_{1} = {F_{1} + F_{2} + F_{3} + F_{4} + F_{5} + F_{6}}} \\ {u_{2} = {F_{2} + F_{3} - F_{5} - F_{6}}} \\ {u_{3} = {F_{1} - F_{4}}} \\ {u_{4} = {{- F_{1}} - F_{3} - F_{5} + F_{2} + F_{4} + F_{6}}} \end{matrix} \right.$

The relationship may be further written as

$\begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \end{bmatrix} = {\begin{bmatrix} 1 & 1 & 1 & 1 & 1 & 1 \\ 0 & 1 & 1 & 0 & {- 1} & {- 1} \\ 1 & 0 & 0 & {- 1} & 0 & 0 \\ {- 1} & 1 & {- 1} & 1 & {- 1} & 1 \end{bmatrix}\begin{bmatrix} F_{1} \\ F_{2} \\ F_{3} \\ F_{4} \\ F_{5} \\ F_{6} \end{bmatrix}}$

Since W is not a square matrix in the above relationship, W⁻¹ may be obtained by means of pseudo-inversion, i.e.

$w^{- 1} = \begin{bmatrix} 0.1667 & 0 & 0.5 & 0 \\ 0.1667 & 0.25 & 0.25 & 0.25 \\ 0.1667 & 0.25 & {- 0.25} & {- 0.25} \\ 0.1667 & 0 & {- 0.5} & 0 \\ 0.1667 & {- 0.25} & {- 0.25} & {- 0.25} \\ 0.1667 & {- 0.25} & 0.25 & 0.25 \end{bmatrix}$

In this way, the power allocation of the six-rotor unmanned aircraft and the corresponding actual motor control quantity may be obtained to control the operation of the motor.

The foregoing description is merely illustrative of preferred embodiments of the present invention, but the scope of protection of the present invention is not limited thereto. Equivalent replacements or modifications made to the inventive concept or technical solution of the present invention by a person skilled in the art within the scope of the disclosure of the present invention fall into the scope of protection of the present invention. 

1. A method for controlling stable flight of an unmanned aircraft, characterized by comprising the steps of: 1) acquiring real-time flight operation data of the aircraft itself, analyzing a kinematic problem of the aircraft, and establishing a dynamics model of the aircraft; 2) constructing a deviation function according to the real-time flight operation data acquired in step 1) and target attitude data, and constructing neurodynamic equations based on the deviation function by using a multi-layer zeroing neurodynamic method, wherein the neurodynamic equations based on the deviation function corresponding to all parameters together constitute a controller of the unmanned aircraft, and output quantities solved from differential equations of the controller are output control quantities of motors of the aircraft; and 3) controlling powers of the motors according to a relationship between the output control quantities solved in step 2) and the powers of the motors of the multi-rotor unmanned aircraft to complete motion control over the unmanned aircraft, specifically: according to a power allocation scheme for the unmanned aircraft, the control quantities solved by the controller have the following relationship with the powers of the motors of the multi-rotor unmanned aircraft: U=WF where U=[u₁ u₂ u₃ u₄]^(T) refers to the output control quantities of the unmanned aircraft, F=[F₁ . . . F_(j)]^(T) refers to the powers of the motors of the unmanned aircraft, j is the number of the motors of the multi-rotor unmanned aircraft, W is a power allocation matrix of the unmanned aircraft, and the matrix W has different forms depending on different structures and the number of rotors, and needs to be determined according to the structure thereof and the number of the rotors; the corresponding powers F of the motors are obtained by means of matrix inversion or pseudo-inversion, i.e.: F=W ⁻¹ U if the matrix W is a square matrix and is reversible, W⁻¹ is obtained by means of an inverse operation, and if W is not a square matrix, W⁻¹ is solved by means of a corresponding pseudo-inverse operation; and the desired powers F of the motors are finally obtained, input voltages of the motors are controlled according to a relationship between the voltages and powers of the motors to control the rotational speeds of the motors, and the control over the powers of the motors is finally realized to complete stable flight control over the unmanned aircraft.
 2. The method for controlling stable flight of an unmanned aircraft according to claim 1, characterized in that the step of analyzing a kinematic problem of the aircraft by a processor mounted thereto specifically comprises: defining a ground coordinate system E and a fuselage coordinate system B, and establishing a relationship E=RB between the ground coordinate system and the fuselage coordinate system by means of a transformation matrix R, where R is expressed as $R = \begin{bmatrix} {\cos\mspace{11mu}\theta\mspace{11mu}\cos\mspace{11mu}\psi} & {\cos\mspace{11mu}\theta\mspace{11mu}\sin\mspace{11mu}\psi} & {{- \sin}\mspace{11mu}\theta} \\ {{\sin\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\cos\mspace{11mu}\psi} - {\cos\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\psi}} & \begin{matrix} {{\sin\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\sin\mspace{11mu}\psi} +} \\ {\cos\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\psi} \end{matrix} & {\sin\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\theta} \\ {{\cos\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\cos\mspace{11mu}\psi} + {\sin\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\psi}} & \begin{matrix} {{\cos\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\sin\mspace{11mu}\psi} -} \\ {\sin\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\psi} \end{matrix} & {\cos\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\theta} \end{bmatrix}$ where ϕ is a roll angle, θ is a pitch angle, and ψ is a yaw angle; ignoring the effect of an air resistance on the aircraft, the stress analysis of the aircraft system in the fuselage coordinate system is as follows ${{\begin{bmatrix} {mI}_{3 \times 3} & 0 \\ 0 & I \end{bmatrix}\begin{bmatrix} \overset{.}{V} \\ \overset{.}{\omega} \end{bmatrix}} + \begin{bmatrix} {\omega \times m\; V} \\ {\omega \times I\;\omega} \end{bmatrix}} = \begin{bmatrix} F \\ \tau \end{bmatrix}$ where m is the total mass of the aircraft, I_(3×3) is a unit matrix, I is an inertia matrix, V is a linear velocity in the fuselage coordinate system, ω is an angular velocity in the fuselage coordinate system, F is a resultant external force, and τ is a resultant torque.
 3. The method for controlling stable flight of an unmanned aircraft according to claim 2, characterized in that the step of establishing a dynamics model of the aircraft specifically comprises: according to the defined ground coordinate system E and fuselage coordinate system B, the relationship E=RB established between the two by means of the transformation matrix R and the stress analysis of the aircraft system in the fuselage coordinate system, obtaining dynamics equations of the multi-rotor aircraft as follows $\left\{ \begin{matrix} {\overset{¨}{x} = \frac{u_{x}u_{1}}{m}} \\ {\overset{¨}{y} = \frac{u_{y}u_{1}}{m}} \\ {\overset{¨}{z} = {{- g} + \frac{\left( {\cos\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\theta} \right)u_{1}}{m}}} \\ {\overset{¨}{\phi} = \frac{{\left( {I_{y} - I_{z}} \right)\overset{.}{\theta}\overset{.}{\psi}} + {lu_{2}}}{I_{x}}} \\ {\overset{¨}{\theta} = \frac{{\left( {I_{z} - I_{x}} \right)\overset{¨}{\psi}\overset{¨}{\phi}} + {lu_{3}}}{I_{y}}} \\ {\overset{¨}{\psi} = \frac{{\left( {I_{x} - I_{y}} \right)\overset{¨}{\phi}\overset{¨}{\theta}} + u_{4}}{I_{z}}} \end{matrix} \right.$ where l is an arm length; g is a gravitational acceleration; x, y, z are respectively position coordinates of the aircraft in the ground coordinate system; {umlaut over (x)}, ÿ, {umlaut over (z)} respectively represent second derivatives of x(t), γ(t) and z(t); ϕ, θ and ψ respectively represent a roll angle, a pitch angle and a yaw angle; {umlaut over (ϕ)}, {umlaut over (θ)} and {umlaut over (ψ)} respectively represent second derivatives of the corresponding parameters; {dot over (ϕ)}, {dot over (θ)}, {dot over (ψ)} respectively represent first derivatives of the corresponding parameters; I_(x), I_(y), I_(z) are respectively rotational inertia of the aircraft in X, Y and Z axes; u_(x)=cos ϕ sin θ cos ψ+sin ϕ sin ψ; u_(y)=cos ϕ sin θ sin ψ−sin ϕ cos ψ; and u₁, u₂, u₃, u₄ are output control quantities.
 4. The method for controlling stable flight of an unmanned aircraft according to claim 1, characterized in that the step of designing a controller of the unmanned aircraft according to a multi-layer zeroing neurodynamic method specifically comprises the steps of: (2-1) designing a deviation function regarding the output control quantity u₁ from the vertical altitude z by means of the multi-layer zeroing neurodynamic method, and designing an altitude controller for the unmanned aircraft according to this deviation function; (2-2) designing a deviation function regarding u_(x) and u_(y) from the horizontal positions x and y by means of the multi-layer zeroing neurodynamic method, designing a position controller for the unmanned aircraft according to this deviation function, and then inversely solving target attitude angles ϕ_(T) and θ_(T); and (2-3) designing a deviation function regarding the output control quantities u₂˜u₄ from the roll angle ϕ, the pitch angle θ and the yaw angle ψ by means of the multi-layer zeroing neurodynamic method, and designing an attitude controller according to this deviation function.
 5. The method for controlling stable flight of an unmanned aircraft according to claim 4, characterized in that the step of designing a deviation function regarding the output control quantity u₁ from the vertical altitude z by means of the multi-layer zeroing neurodynamic method, and designing an altitude controller for the unmanned aircraft according to this deviation function specifically comprises: for the vertical altitude z, according to the target altitude value z_(T) and the actual altitude value z in the Z axis, defining a deviation function as E _(z) ={umlaut over (z)}−{umlaut over (z)} _(T)+2γ(ż−ż _(T))+γ²(z−z _(T))  (9) according to the dynamics equations of the aircraft, simplifying (9) into E _(z) =a _(z) u ₁ +b _(z)  (10) where ${a_{z} = \frac{\cos\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\theta}{m}},{b_{z} = {{- g} - {\overset{¨}{z}}_{T} + {2{\gamma\left( {\overset{.}{z} - {\overset{.}{z}}_{T}} \right)}} + {\gamma^{2}\left( {z - z_{T}} \right)}}},$ and γ is a constant; and obtaining its derivative as follows Ė _(z) =a _(z) {dot over (u)} ₁ +{dot over (a)} _(z) u ₁ +{dot over (b)} _(z)  (11) using the multi-layer zeroing neurodynamic method to design Ė _(z) =−γE _(z)  (12) substituting equations (10) and (11) into equation (12) and perform collating to obtain a _(z) {dot over (u)} ₁=γ(a _(z) u ₁ +b _(z))−{dot over (b)} _(z) −{dot over (a)} _(z) u ₁  (13).
 6. The method for controlling stable flight of an unmanned aircraft according to claim 4, characterized in that the step of designing a deviation function regarding u_(x) and u_(y) and a position controller for the unmanned aircraft specifically comprises: for the horizontal position x, according to the target value x_(T) and the actual value x in the X axis, defining a deviation function as E _(x) ={umlaut over (x)}−{umlaut over (x)} _(T)+2γ({dot over (x)}−{dot over (x)} _(T))+γ²(x−x _(T))  (22) according to the dynamics equations of the aircraft, simplifying equation (22) into E _(x) =a _(x) u _(x) +b _(x)  (23) where ${a_{x} = \frac{u_{1}}{m}},$ and b_(x)=−{umlaut over (x)}_(T)+2γ({dot over (x)}−{dot over (x)}_(T))+γ²(x−x_(T)); and obtaining its derivative as follows Ė _(x) =a _(x) {dot over (u)} _(x) +{dot over (a)} _(x) u _(x) +{dot over (b)} _(x)  (24) using the multi-layer zeroing neurodynamic method to design Ė _(x) =−γE _(x)  (25) substituting equations (23) and (24) into equation (25) and perform collating to obtain a _(x) {dot over (u)} _(x)=−γ(a _(x) u _(x) +b _(x))−{dot over (b)} _(x) −{dot over (a)} _(x) u _(x)  (26) for the horizontal position y, according to the target value y_(T) and the actual value y in the Y axis, defining a deviation function as E _(y) =ÿ−ÿ _(T)+2γ({dot over (y)}−{dot over (y)} _(T))+γ²(y−y _(T))  (35) according to the dynamics equations of the aircraft, simplifying equation (35) into E _(y) =a _(y) u _(y) +b _(y)  (36) where ${a_{y} = \frac{u_{1}}{m}},$ and b_(y)=−ÿ_(T)+2γ({dot over (y)}−{dot over (y)}_(T))+γ² (y−y_(T)); and obtaining its derivative as follows Ė _(y) =a _(y) {dot over (u)} _(y) +{dot over (a)} _(y) u _(y) +{dot over (b)} _(y)  (37) using the multi-layer zeroing neurodynamic method to design Ė _(y) =−γE _(y)  (38) substituting equations (36) and (37) into equation (38) and perform collating to obtain a _(y) {dot over (u)} _(y)=−γ(a _(y) u _(y) +b _(y))−{dot over (b)} _(y) −{dot over (a)} _(y) u _(y)  (39).
 7. The method for controlling stable flight of an unmanned aircraft according to claim 6, characterized in that the calculation formulas of inversely solving the target attitude angles ϕ_(T) and θ_(T) are: u_(x) and u_(y) solved from position controller equations (26) and (39) are $\left\{ \begin{matrix} {u_{x} = {{\cos\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\cos\mspace{11mu}\psi} + {\sin\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\psi}}} \\ {u_{y} = {{\cos\mspace{11mu}\phi\mspace{11mu}\sin\mspace{11mu}\theta\mspace{11mu}\sin\mspace{11mu}\psi} - {\sin\mspace{11mu}\phi\mspace{11mu}\cos\mspace{11mu}\psi}}} \end{matrix} \right.$ so that the inversely solved target angle values ϕ_(T) and θ_(T) are $\begin{matrix} {\left\{ \begin{matrix} {\phi_{T} = {\sin^{- 1}\left( {{u_{x}\mspace{11mu}\sin\mspace{11mu}\psi} - {u_{y}\mspace{14mu}\cos\mspace{11mu}\psi}} \right)}} \\ {\theta_{T} = {\sin^{- 1}\left( \frac{u_{x} - {\sin\mspace{11mu}\phi_{T}\mspace{14mu}\sin\mspace{11mu}\psi}}{\cos\mspace{11mu}\phi_{T}\mspace{14mu}\cos\mspace{11mu}\psi} \right)}} \end{matrix} \right..} & (40) \end{matrix}$
 8. The method for controlling stable flight of an unmanned aircraft according to claim 4, characterized in that the step of designing a deviation function regarding the output control quantities u₂˜u₄ from the roll angle ϕ, the pitch angle θ and the yaw angle ψ by means of the multi-layer zeroing neurodynamic method, and designing an attitude controller according to this deviation function specifically comprises: for the roll angle ϕ, according to the target angle ϕ_(T) solved in (40) and the actual angle ϕ, defining a deviation function as E _(ϕ)={umlaut over (ϕ)}−{umlaut over (ϕ)}_(T)+2γ({dot over (ϕ)}−{dot over (ϕ)}_(T))+γ²(ϕ−ϕ_(T))  (49) according to the dynamics equations of the aircraft, simplifying equation (49) into E _(ϕ) =a _(ϕ) u ₂ +b _(ϕ)  (50) where ${a_{\phi} = \frac{l}{I_{x}}},{{{{and}\mspace{14mu} b_{\phi}} = {\frac{\left( {I_{y} - I_{z}} \right)\overset{.}{\theta}\overset{.}{\psi}}{I_{x}} - {\overset{¨}{\phi}}_{T} + {2{\gamma\left( {\overset{.}{\phi} - {\overset{.}{\phi}}_{T}} \right)}} + {\gamma^{2}\left( {\phi - \phi_{T}} \right)}}};}$ and obtaining its derivative as follows Ė _(ϕ) =a _(ϕ) {dot over (u)} ₂ +{dot over (a)} _(ϕ) u ₂ +{dot over (b)} _(ϕ)  (51) according to the multi-layer zeroing neurodynamic method, designing Ė _(ϕ) =−γE _(ϕ)  (52) substituting equations (50) and (51) into equation (52) and perform collating to obtain a _(ϕ) {dot over (u)} ₂=−γ(a _(ϕ) u ₂ +b _(ϕ))−{dot over (b)} _(ϕ) −{dot over (a)} _(ϕ) u ₂  (53) for the pitch angle θ, according to the target angle θ_(T) solved in (40) and the actual angle θ, defining a deviation function as E _(θ)={umlaut over (θ)}−{umlaut over (θ)}_(T)+2γ({dot over (θ)}−{dot over (θ)}_(T))+γ²(θ−θ_(T))  (62) according to the dynamics equations of the aircraft, simplifying equation (62) into E _(θ) =a _(θ) u ₃ +b _(θ)  (63) where ${a_{\theta} = \frac{l}{I_{y}}},{{{{and}\mspace{14mu} b_{\theta}} = {\frac{\left( {I_{z} - I_{x}} \right)\overset{.}{\phi}\overset{.}{\psi}}{I_{y}} - {\overset{¨}{\theta}}_{T} + {2{\gamma\left( {\overset{.}{\theta} - {\overset{.}{\theta}}_{T}} \right)}} + {\gamma^{2}\left( {\theta - \theta_{T}} \right)}}};}$ and obtaining its derivative as follows Ė _(θ) =a _(θ) {dot over (u)} ₃ +{dot over (a)} _(θ) u ₃ +{dot over (b)} _(θ)  (64) according to the multi-layer zeroing neurodynamic method, designing Ė _(θ) =−γE _(θ)  (65) substituting equations (63) and (64) into equation (65) and perform collating to obtain a _(θ) {dot over (u)} ₃=−γ(a _(θ) u ₃ +b _(θ))−{dot over (b)} _(θ) −{dot over (a)} _(θ) u ₃  (66) for the yaw angle ψ, according to an artificially set angle ψ_(T) and the actual angle ψ, defining a deviation function as E _(ψ)={umlaut over (ψ)}−{umlaut over (ψ)}_(T)+2γ({dot over (ψ)}−{dot over (ψ)}_(T))+γ²(ψ−ψ_(T))  (75) according to the dynamics equations of the aircraft, simplifying equation (75) into E _(ψ) =a _(ψ) u ₄ +b _(ψ)  (76) where ${a_{\psi} = \frac{1}{I_{z}}},{{{{and}\mspace{14mu} b_{\psi}} = {\frac{\left( {I_{x} - I_{y}} \right)\overset{.}{\phi}\overset{.}{\theta}}{I_{z}} - {\overset{¨}{\psi}}_{T} + {2{\gamma\left( {\overset{.}{\psi} - {\overset{.}{\psi}}_{T}} \right)}} + {\gamma^{2}\left( {\psi - \psi_{T}} \right)}}};}$ and obtaining its derivative as follows Ė _(ψ) =a _(ψ) {dot over (u)} ₄ +{dot over (a)} _(ψ) u ₄ +{dot over (b)} _(ψ)  (77) according to the multi-layer zeroing neurodynamic method, designing Ė _(ψ) =−γE _(ψ)  (78) substituting equations (76) and (77) into equation (78) and perform collating to obtain a _(ψ) {dot over (u)} ₄=γ(a _(ψ) u ₄ +b _(ψ))−{dot over (b)} _(ψ) −{dot over (a)} _(ψ) u ₄  (79).
 9. The method for controlling stable flight of an unmanned aircraft according to claim 4, characterized in that the step in which the designed altitude controller, position controller and attitude controller together constitute a stable aircraft of the multi-rotor unmanned aircraft specifically comprises: a _(z) {dot over (u)} ₁=−γ(a _(z) u ₁ +b _(z))−{dot over (b)} _(z) −{dot over (a)} _(z) u ₁  (13) a _(ϕ) {dot over (u)} ₂=−γ(a _(ϕ) u ₂ +b _(ϕ))−{dot over (b)} _(ϕ) −{dot over (a)} _(ϕ) u ₂  (53) a _(θ) {dot over (u)} ₃=−γ(a _(θ) u ₃ +b _(θ))−{dot over (b)} _(θ) −{dot over (a)} _(θ) u ₃  (66) a _(ψ) {dot over (u)} ₄=−γ(a _(ψ) u ₄ +b _(γ))−{dot over (b)} _(ψ) {dot over (a)} _(ψ) u ₄  (79) obtaining a controller of the unmanned aircraft according to equations (13), (53), (66) and (79), wherein the controller can be implemented by a network structure; the controller of the unmanned aircraft is capable of controlling the stable flight of the unmanned aircraft; and the controller is written in the following form: $\left\{ \begin{matrix} {{\overset{.}{u}}_{1} = \frac{{- {\gamma\left( {{a_{z}u_{1}} + b_{z}} \right)}} - {\overset{.}{b}}_{z} - {{\overset{.}{a}}_{z}u_{1}}}{a_{z}}} \\ {{\overset{.}{u}}_{2} = \frac{{- {\gamma\left( {{a_{\phi}u_{2}} + b_{\phi}} \right)}} - {\overset{.}{b}}_{\phi} - {{\overset{.}{a}}_{\phi}u_{2}}}{a_{\phi}}} \\ {{\overset{.}{u}}_{3} = \frac{{- {\gamma\left( {{a_{\theta}u_{3}} + b_{\theta}} \right)}} - {\overset{.}{b}}_{\theta} - {{\overset{.}{a}}_{\theta}u_{3}}}{a_{\theta}}} \\ {{\overset{.}{u}}_{4} = \frac{{- {\gamma\left( {{a_{\psi}u_{4}} + b_{\psi}} \right)}} - {\overset{.}{b}}_{\psi} - {{\overset{.}{a}}_{\psi}u_{4}}}{a_{\psi}}} \end{matrix} \right.$ a multi-layer zeroing neural network is constructed from the differential equations of the controller, and the control quantities of the unmanned aircraft are solved by means of the multi-layer zeroing neural network. 